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Abstract 

We show that the equation in the title (with $„ the nth cyclotomic 
polynomial) has no integer solution with n > 1 in the cases (m,p) = 
(15, 41), (15, 5581), (10, 271). These equations arise in a recent group 
theoretical investigation by Z. Akhlaghi, M. Khatami and B. Khosravi. 

1 Introduction 

In the recent work [TJ by Zeinab Akhlaghi, Maryam Khatami and Behrooz 
Khosravi, some Diophantine equations come up in a group theoretical con- 
text. In particular, Zeinab Akhlaghi posed the following problems to us. 

• Which primes P of the form P = 2 ■ Al 2a — 1 can also be written as 
P = ^15 (f?) , with q a prime power? 

• Which primes P of the form P = 2 • 5581 2a — 1 can also be written as 
P = <&is(±q), with q a prime power? 

• Which primes P of the form P = 2 ■ 271 2a — 1 can also be written as 
P = $io(g 2 ), with q a prime power? 

Here $ m is the mth cyclotomic polynomial. In particular, 

$ 15 (x) = x 8 - x 7 + x 5 - x 4 + x 3 - X + 1, 
$io(x) = x A - x 3 + x 2 - X + 1. 

Note that <S> w (q 2 ) = <S>2o(q)- 
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Typical for Diophantine equations arising in group theory is the oc- 
currence of primes, and so the above present some 'typical' examples of 
equations so arising. 

Given a group G, let 7r(G) denote the set of primes q such that G contains 
an element of order q. Then the prime graph r(G) of G is defined as the 
graph G with vertex set 7r(G) in which two distinct primes q, q' £ 7r(G) are 
adjacent if G contains an element of order qq'. Akhlaghi et al. [1] show, 
using Theorem Q] below and various already known Diophantine results, 
that in case p is an odd prime and k > 1 is odd, then PGL(2,p fc ) is uniqely 
characterized by its prime graph, i.e. there is no other group having the 
same prime graph. 

In this paper we will prove the following result, implying the answer 
"None" to the above three problems, 

Theorem 1. Let (m,p) = (15, 41), (15, 5581), or (10,271). Then the Dio- 
phantine equation 

3> m {x) + l = 2p n (1) 
has no integer solutions (n,x) with n>l. 

In the literature, by different methods, some equations of a similar nature 
have been studied, e.g. the equations <&m{x) = P n and <fr m (x) = p n + 1, with 
m a prime. The first equation is a special case of the Nagell-Ljunggren 
equation (x m — l)/(x — 1) = p n and is studied in many papers (for a survey 
see [2]). For a non-existence result of solutions of the second equation see 
Lei]. 

General results on solutions of equations of the form fix) = by m (see the 
book by Shorey and Tijdeman [8J) imply that for an arbitrary, but fixed m > 
3, equation ([I]) has finitely many solutions (x,p,n) with max{|x|,p, re} < C, 
with C a computable number. Formulated in this generality, C which comes 
from applying the theory of linear forms in logarithms, will be huge. 

In Section [21 we give an elementary proof of a lower bound n > 239, and 
a related heuristic argument why we do not expect any solutions for these 
problems. In Section El we use algebraic number theory and a deep result 
from transcendence theory to deduce an upper bound n < 2.163 • 10 27 for 
re satisfying (pQ) . Then in Section [5] the LLL algorithm will be invoked to 
efficiently reduce this bound to n < 59. In this way we obtain a rigorous, 
albeit computational, proof of Theorem [TJ We note that our method should 
work in principle for other equations of the type f{x) = ap n , when / is a 
fixed polynomial with integral coefficients and at least three distinct roots, 
a > 1 is a fixed integer, and p is a fixed prime not dividing the discriminant 
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of /. The nature of our method is algorithmic in the sense that for every 
single choice of parameters the details of the method have to be worked 
through separately. 

An extended version of this paper including an appendix with numerical 
material related to the application of the LLL-method and a subsection on 
reducing the equation modulo a prime q ^ p (a subsection not relevant for 
the proof of our main result), is available as MPIM preprint 2009-62 [5|. 

2 j9-adic considerations 

Without loss of generality we may assume that \x\ > 2 and n > 1. We write 
fm{x) = &m(x) + 1 and d = deg/ m for m = 10,15. Elementary calculus 
shows that for all x 

(\ X \ - if < $ m (x) < f m (x) < {\x\ + l) d . (2) 

See e.g. [3] for some similar estimates. We start with seeing what information 
we can derive from studying the p-adic roots of f m . If (x, n) is a solution of 
(PQ), then there is a root 

oo 

x = Yj a kP k (with dfc G {0, 1, . . . ,p — 1}) 

k=0 

of f m in Qp such that x = x (mod p n ). Note that if ao 7^ then the p-adic 
expansion of — x is 

00 

-x = (p - a ) + - 1 - a k )p k . 

k=l 

Now (j2J) with f m (x) = 2p n implies that 

\x\ < 2 l ' d p n ' d + l< 2p n l d , 

and this immediately implies that, in the case x > 

a k = for all k G N with [(n + l)/dj + 1 < k < n - 1, 

and in the case x < 

a fc = p - 1 for all k G N with [(n + l)/d] + 1 < fc < n - 1. 

In other words, the existence of a solution n of ([T]) implies that of the first n 
p-adic digits of the root x, the last consecutive ~ n(l — 1/d) all have to be 
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equal to or p — 1, in respectively the cases x > and x < 0. This seems 
unlikely to happen, as can easily be verified experimentally for not too large 
n. It seems not unreasonable to expect that the p-adic digits of the roots x 
are uniformly distributed over {0, 1, ... ,p — 1}, and that these distributions 
per digit are independent. Then the probability that n(l — 1/d) specific 
consecutive digits are all respectively p — 1 is p~ nt ~ l ~ l l d ) ; an d the expected 
number of solutions is at most 



We conclude that if p is large, then likely there are no solutions. If p is small 
and there is no solution with n small, then very likely there are no solutions 
at all. 

A minor variation of the above argument suggests that in case d > 3 
there are only finitely many solutions (x,p,n) of (H|) with n > 2. As we 
already remarked in the introduction, this result is known to be true, see 



Explicit computation of the »-adic root x up to some finite precision is a 
quick way to rule out small values of n. We now give details for the cases that 
are of interest to us. Note that p-adic roots of polynomials are quite easy 
to compute by Hensel lifting (i.e. the p-adic version of the Newton- Raphson 
method). 

In the case p = 41 there is one 41-adic root of /15. Its sequence of 41-adic 
digits is 

8, 18, 3, 17, 9, 14, 12, 38, 31, 35, 19, 25, 19, 38, 25, 24, 1, 18, 
25, 10, 14, 29, 31, 18, 36, 2, 24,... 

The smallest k such that = or 40 is k = 53. Hence a solution of ([I]) 
implies k > 53, which in turn implies \_(n + l)/8j +1 > 53, so n > 415. 

Two remarks are in place. Firstly, we did not even bother to use con- 
secutive zeros, we used only one. Indeed, 054 = 15, so we could sharpen 
our result easily. But we have to stop somewhere, and the result n > 415 is 
sufficient for the moment. And secondly, it should be noted that the com- 
plexity of this method is exponential, as to compute the nth p-adic digit 
we have to compute with numbers of the size p n . This makes this method 
unrealistic for values of n that become larger than a few thousand. 

In the case p = 5581 there are two 5581-adic roots of /15. Their sequences 



n=l 
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of 5581-adic digits are 



257, 64, 5438, 1453, 629, 833, 3090, 5096, 4809, 1493, 4462, 1922, 
4807, 782, 3819, 2190, 99, 2554, 3603, 4471, 1034, 1407, 3688, ... 

and 

4477, 3993, 3590, 3157, 3667, 3404, 2233, 3440, 3784, 2333, 900, 
2522, 184, 1707, 5103, 2005, 5325, 1780, 4765, 2645, 3577, ... 

In both cases we computed up to k = 502 and did not encounter a or a 
5580. As above it follows that n > 4015. 

In the case p = 271 there is one 271-adic root of /io- Its sequence of 
271-adic digits is 

241, 8, 147, 250, 135, 263, 1, 126, 89, 262, 149, 20, 147, 78, 
220, 219, 176, 148, 206, 255, 38, 115, 186, 178, 235, ... 

The smallest k such that a& = or = 270 is k = 61. Hence a solution of 
([1]) implies k > 61, which in turn implies [(n + 1)/4J + 1 > 61, so n > 239. 

Using the above results we infer that on heuristic grounds with prob- 
ability at most io -1000 equation ([1]) has a non-trivial solution. Since in 
mathematics one has to prove assertions beyond 'unreasonable doubt', we 
cannont conclude our paper at this point. 

3 Finding an upper bound 

We start with giving some data on relevant algebraic number fields. Then 
we derive from equation ([1]) an S'-unit inequality, to which we apply tran- 
scendence theory to find an explicit upper bound for n. 



3.1 Field data 



We have 



/is(x) 
fio(x) 



x 8 - x 7 + x 5 - x 4 + x 3 - x + 2, 
x 4 - x 3 + x 2 - x + 2, 



and we write equation (JXJ) as 



f m (x) = 2p n , 



(3) 
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where (m,p) = (15, 41), (15, 5581), (10, 271). For brevity we will refer to 
these cases as the cases p = 41,5581,271 or m = 15,10. In Section [21 we 
have seen that n > 239, and we may assume that \x\ > 2. 

The polynomials f m are irreducible and have no real roots. We label the 
roots as follows: 



root 


fib 


ho 


qW 


1.0757 


. .+0.4498 


. . . i 


0.9734... +0.7873... i 


a( 2 ) 


0.6243 


. .+0.8958 


. . . i 


-0.4734... +1.0255... i 


a( 3 ) 


-0.1701 


..+1.0292 


. . . i 






-1.0299 


..+0.2698 


. . . i 






a^ = aSi~ 


-4) f or j — 


5,6,7,8 


a (j) = for j _ 3j 4 


max \a*j) < 


1.167 


1.252 



We write K m for the field Q(a) where a is a root of f m {x) = 0, so that 
d = deg f m = [K m : Q], i.e. d = 8 for m = 15 and d = 4 for m = 10. 

We need a lot of data on these fields. We used Pari [7] to obtain the 
data given below. 

The discriminants of Ki 5 ,Ki are respectively 682862912 = 2 6 -83-128551 
and 1396 = 2 2 • 349. In both cases a generates a power integral basis. 
Fundamental units are: 

for m = 15 : 0i = a 7 + a 4 + a 2 + a — 1, 
f3 2 = a 6 - a 5 + a 4 + a - 1, 
/?3 = a 2 — a + 1, 

for m = 10 : fix = a 3 - a 2 + 1. 

The regulators are 4.2219 1.1840 .. . respectively. The class groups of 
both fields is trivial. 

The prime decomposition of 2 is 

for m = 15 : 2 = a(a + l) 4 (a 3 - a 2 + l)/?r 2 /? 2 , 
form =10: 2 = -a(a - l) 3 /3f 1 . 

Thus the prime ideals of norm 2 are (a), (a+1) when m = 15, and (a), (a— 1) 
when m = 10. 

The prime decomposition of p in the field K m is as follows: 
for p = 41: 41 = 7172, where 

71 = - a 7 + a 6 + a 5 - 2a 4 + a 3 + a 2 - a + 1, ^(71) = 41, 

72 = _4a 7 + 13a 6 - 19a 5 + 8a 4 - 14a 3 + 7a 2 + 15a + 1, N{j 2 ) = 41 7 , 
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for p = 5581: 5581 = 71727374, where 

71 = a 6 -a 5 -2a + l, N(m) 

72 = 2a 5 + a 2 + a + 1, iV(7 2 ) 

73 = -3a 7 - a 6 + 7a 5 - 4a 4 - 5a 3 + 7a 2 + a + 1, iV( 73 ) 

74 = 85a 7 - 41a 6 - 112a 5 + 55a 4 - 21a 3 + 

+134a 2 + 92a - 135, #(74) 

for p = 271: 271 = 7172, where 

71 = -2a 3 + 4a 2 - 4a + 3, iV( 7 i) = 271, 

72 = -18a 3 + 16a 2 + 44a + 53, JV(^j) = 271 3 . 

3.2 Deriving an S*-unit inequality 

If x is an integer satisfying ([3|), then it follows that in Ojj we have 

(x - a)z = 2p n 

for a z S Ok- Thus, we can write (taking 73 = 74 = 1 in the cases p = 
41,271) 

x -a = «57r i 7 2 n2 7 3 n3 7 4 n4 /?, * = (2/<5) 7 r ni 7 2 n ^it^ l2~ ni 

where 5 | 2 and is a unit. Taking norms we find 

2p n = N{x - a) = N(5)p Cini+C2n2+C3ns+C4n \ 

where (01,02,03,04) = (1, 7, 0, 0), (1, 1, 2, 4), (1, 3, 0, 0) for respectively p = 
41, 5581, 271. It follows that N(5) = 2, and n = cini + c 2 n-2 + 037x3 + 04/14. 

First observe that < nj < n is impossible. Indeed, for if not, then there 
exists k S {1,2,3,4} such that 7^ 7^ 1 divides both x — a and z. Observe 
that if a = a'*', then z = Ylj^x — a^). Thus, if p is some prime ideal of 
Ojg- dividing 7^, then p divides both x — a' 1 ' and x — a^> for some j ^ i. 
In particular, p divides and thus also A(/ m ). Since this last 

number is an integer and p has norm a power of p in K m , it would follow 
that p divides A(/ m ), which is not the case. Thus, the only possibilities 
are n« € {0, n} for all i. The equation n = cini + C2n 2 + 03723 + 0477,4 
now has only the solutions (t7i,77 2 ) = (n, 0) in the cases p = 41,271, and 
(774, ?i 2 , 713, 714) = (n, 0, 0, 0), (0, 77, 0, 0) in the case p = 5581. 

We get the following equations: 

p = 41 : x - a = ±5~f n (3™ 1 (3™ 2 (3™ 3 , 5 = a,a + l, 7 = 71, 

p = 5581: x-a = ±8-f n p™ 1 p% l2 P™ 3 , 6 = a,a + l, 7 = 7^72, (4) 

p = 271: x - a = ±5-y n (3™ 1 , <5 = a,a-l, 7 = 71. 



= 5581, 
= 5581, 
= 5581 2 , 

= 5581 4 , 
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Now we could proceed by conjugating equation and eliminating x to get 
a unit equation. However, this resulting unit equation will live in the field 
Q[a, a], which is of degree d(d — 1), because the Galois group of f m (x) over 
Q is Sd- Since estimates for linear forms in logarithms are quite sensitive 
to the degree, we will continue to work in K m . We proceed as follows. For 
convenience in the cases p = 41, 271 we put fa = fa = 1 and ni2 = = 0. 
We have from (JI|) that 



2p n 



x — a 



(5) 



Putting y = x — a, Taylor's formula yields z = > - — hence 



i=l 



z-y 



1-1 



ST^ Jm («) j-1 



1=1 



Let us now make some estimates. Observe that the lower bound n > 239 
from Section [2] is amply sufficient to guarantee p n > (2 • 10 10 ) d . Then ([2]) 
implies 



\y\ 



\x-a\ > \x\-\a\ > f m (xfl d -l-\a\ > 2 1 ' d p n ' d -2.2h2 > C lP n / d , (6) 



where C\ = 1.090 when m = 15 and C\ = 1.189 when m = 10. Hence, 



y\ > 2 ■ 10 10 . We now compute upper bounds for 



getting 





1 


2 


3 4 5 6 7 


l/i ( 5 j («)IA! < 


16.40 
6.977 


56.37 
9.261 


109.6 126.7 90.07 39.00 9.489 
5.021 



so that 



z-y 



d-l 



< \y 



id— 2 



d-l 

£ 

i=l 



/(i) 



d-2 



where C2 = 9.490 for m = 15 and C2 = 5.022 for m = 10, because \y\ > 
2 ■ 10 10 . Thus, 

(7) 



1 



If 



d-l 



c 2 c 3 
|y| p d 
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c 2 

where C3 > — , so C3 = 8.706 for m = 15 and C3 = 4.223 for m = 
Using equations (j4|), ([5]) and ([7]), we get the S"-unit inequality we want: 



1 



P 



T 



— 8mi fl-8m2 0—81113 
PS 



< 



c 3 

n/d' 



P 



10. 



(8) 



3.3 Applying transcendence theory 

We shall apply a linear form in logarithms to bound the expression on the 
left of inequality © from below. We first check that it is not zero. If it 
were, then since it comes from rewriting the left hand side of inequality ([7]), 
we would get that z = y^ 1 ■ Since yz = 2p n , we get that y d = 2p n , which 
violates the prime decomposition of 2 in K m . 

Next, we need to bound m\,vri2 and 711,3 m terms of n. Since p n l d > 
2 • 10 10 , it follows from that 

\ l / d 4- 1 4- Irvl <- 9 1 M n ™/ c( 4-9 9.^9 S r.,-n n l d 



\v\ 



\x— a\ < \x\ + \a\ < f n 



XH-+1+H < 2 i /> n/a +2.252 < C 4 p n/a , (9) 



where C4 = 1.091 for m = 15 and C4 = 1.190 for m = 10. Now we take 
absolute values of the conjugates of equation @, and rewrite them as 



\x — aW 



7i 



(0 



(i) 



« 



(0 



(10) 



We computed: 



forp = 41: 0.2714 < |<yW I < 2.124, 0.5676 < |-y^)| < 5.349, 
for p = 5581 : 0.2714 < jtfW | < 2.124, 1.522 < I7W | < 5.531, 
for p = 271: 0.7877 < |«jW | < 1.796, 2.253 < j-y^j < 7.307, 



and thus 




pl/d 


,log 


max 7W 


min 7M 


pl/d 


c 4 


log 


max |<5^ 


min |5W 


Ci 




where for p = 41 we have C5 = 1.213, Cg = 1.392, for p = 5581 we have 
C 5 = 0.6584, C 6 = 1.392, and for p = 271 we have C 5 = 0.5884, C 6 = 0.4126. 
It follows from © and © that 



log 



(*) 



\x — a 



< C 5 n + C 6 . 
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Writing ui for the logarithm of the left hand side of equations (fT0|) , we get 
that 



Ui = m\ log 



(0 



m 2 log 



+ m 3 log 



for three conjugates i, (11) 



?(0 



and hence < C^n + C% for all i. If m = 10 this simply states log j3\ 

1.184 (this is the regulator of K\q), as then 02 = 03 = 1, and thus \m 
{C§n + Ce)/1.184. If m = 15, solving the system (jlip with Cramer's rule, 
we get that 

3(C 5 n + C 6 )R 2 



max{|mi|, \m2\, \ms \ } < 



Rf: 



where R 2 is the maximal absolute value of all the 2x2 minors of the co- 
efficient matrix appearing in formula (jlip whose determinant is Rg. The 
minor largest in absolute value is the (2, 1) minor obtained by eliminating 
the second row and first column, and its value is R 2 < 2.746. Putting all 
this together gives 

max{|mi|,|m 2 |,|m 3 |}<C 7 n + C8, 

where C 7 = 2.369, C 8 = 2.718 when p = 41, C 7 = 1.286, C 8 = 2.718 when 
p = 5581, and C 7 = 0.4970, C 8 = 0.3485 when p = 271. 

The next step is to prepare for the application of a deep result from 
transcendence theory. We return to inequality (JSj) and rewrite it as 



i=l 



< 



p 



n/d' 



(12) 



where r = 5 when m = 15 and r = 3 if m = 10, and 



and b\ 
ing 



m 



i,6 2 



5 d> 



V2 



= 2. 

-dmi, 64 



T]3 = 01, m = 02, 115=03, 



-dm,2,b5 = —dm^ are integers satisfy- 



B = max \bi\ < d(C 7 n + Cg). 
Recall that for an algebraic number ij having 

d 

aol[(X-r,®) 

i=l 
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as minimal polynomial over the integers, the logarithmic height is defined 
as 

1 / d 

h(rj) = — ( log | a-o | + log max 

\ i=l 

With this notation, Matveev [6] proved the following deep theorem. 

Theorem 2. Let IK be a field of degree D, 771, .. . ,rjk be nonzero elements 
of K, and b\, . . . , integers. Put 

B = max{|6i|, . . . , \b k \} 

and 

k 

A=i-rR- 

1=1 

Let Ax , . . . , Af; be real numbers such that 

Aj > m&x{Dh(r]j), \ log rjj |, 0.16}, j = 1, . . . ,k. 
Then, assuming that A / 0, we have 

k 

log |A| > -3 • 30 k+4 (k + lf 5 D 2 {l + log D)(l + log{kB)) A, t . 

i=i 

We apply Matveev's result to get a lower bound on the expression ap- 
pearing in the left hand side of (|12p with k = r + 1. We take the field to be 
our K m , so D = d. We also take 7^, bi as in (fT2|) . 

We computed as leading coefficients oq of minimal polynomials: 



m 


5 


m 




m 






15 


a 


a 


= 2< 


a 


= P 7 


a = 1 




a + 1 


a 


= 2 4 








10 


a 


a 


= 2 3 


a 


= p A 


a = 1 




a — 1 


a 


= 2 









and for the Aj we found 



p 


A, < 


A 2 < 


^3 < 


A 4 < 


A 5 < 


41 


25.02 


47.80 


4.371 


4.247 


2.976 


5581 


25.02 


74.22 


4.371 


4.247 


2.976 


271 


3.988 


21.52 


2.634 







,(0 



> 
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Thus, by Matveev's bound we have that 

|logA| > -Cg(l + log(rS)), 
where C 9 > 3 • 30 r+4 (r + l) 5 - 5 d 2 (l + logd)AiA 2 . . . A r satisfies 



for p = 41 
for p = 5581 
for p = 271 



C 9 = 1.465 ■ 10 25 , 
C 9 = 2.275 • 10 25 , 
C 9 = 1.160 • 10 18 . 



Comparing this with the fact that B < d(C 7 n + Cg) and with inequality 
(USD, we get 



logp 



Concretely: 
for p = 41 

for p = 5581 

for p = 271 



-n - log C 3 < - log |A| < C 9 (l + \og{rd(C 7 n + C 8 ))). 

0.4641n - 2.165 < 1.465 • 10 25 (1 + log(94.80n + 108.8)) 
implying n < N = 2.163 • 10 27 , 

1.078n - 2.165 < 2.275 • 10 25 (1 + log(51.45n + 108.8)) 
implying n < N = 1.424 • 10 27 , 

1.400n - 1.441 < 1.160 • 10 18 (1 + log(5.964n + 4.182)) 
implying n< N = 3.970 • 10 19 . 



4 Reducing the upper bound 

So, it remains to solve 

r 

-drrii-2 



1 



n 

i=3 



V; 



max m, 



< 



C3_ 

n n/d'' 



< d(C 7 n + C 8 ), 
n < N. 

This is a finite problem, but the upper bound N is way too large to apply 
brute force or the method from Section [2j Efficient methods for solving such 
problems based on lattice basis reduction using the LLL algorithm exist, see 
[9], and they work quite well in our case. Here are the details. 
We put 



A ' 



0) 



J log 




for 


% = 1,2 




\ -d\og 


4 f) 


for 


i = 3, ... ,r. 


3 



1, 



1. 
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Let 



00 



mi A3 + . . . + m r _2A^ for j = 1, . . . , r — 1. 



By ()12p . the real linear forms A"'' satisfy 



A (i) 



< 



1 - e 



A (i) 



< 



n 

i=3 



,0') 



-drrii- 



< 



n/d' 



(13) 



We let if be some constant slightly larger than N^~ l ''^ r ~ 2 \ i.e. iV 4//3 when 

m = 15 and r = 5, and iV 2 when m = 10 and r = 3. We write #^ = 

for i = 1, . . . , r, where [•] denotes rounding to the nearest integer. We put 

(A') (i) = e[ j) + n9^ j) + miflf + . . . + m r _ 2 ^ j) - 

Then 

ifA (j) - (A') (j) 



1 n r — 2 
< - + - H — max|mi| < Ci n + Cn, 



where Cio = \ + L ^dCj and Cn = ^ + ^-^dC^. Then n > N implies 



(A') (i) 



A 0) 



+ c 10 n + c 



11- 



(14) 



We now look at the matrix T and the vector y given as 



for m = 15 



#0') gu) gu) 
.(*) fl (fc) fl (*) 

V o o o i / 



V o J 



where k) € {(1, 2, 3), (1, 2, 4), (1, 3, 4), (2, 3, 4)}, 



for m = 10 



1 



3 « N 



where i € {1, 2}. 
Observe that for x = (mi, . . . , ?n r _2> n ) T 

rx-y=((A')«,(A') W ,(A') (fc) ,n)' 



resp. 



(A')« 



7? 
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The columns of T generate a sublattice of Z r 2 . Let d(T, y) = min \Tx — y\ 
be the distance from y to the nearest lattice point. From (fT4"|) we find 



d(T,y) < \Tx-y\ < J (r - 2) (K max | \U) | + C 10 N + C u ) 2 + N 2 . (15) 



Put 



ArV(r-2) / U(T,y) 2 -N 2 



iiT V r-2 



If c happens to be a positive real number, then combining ()13p and (I15p we 
get for A = A^), such that |A| = max |A^)| satisfies 

ciV -V(r-2) < | A | < ^3 



pn/d ' 

and hence 

n < I log C 3 — log c H log N 

logp \ r — 2 

In particular, if c is reasonable, that is, not too tiny, then the above bound is 
a reduced upper bound for n. We can argue that this is reasonable, because 
if the lattice is generic, that is, if it satisfies 

d(T,y) « det(r) 1 / dimr « K (r- 2 )/(r-D ^ 

then with the choice of K being somewhat larger than jv( r ~ 1 )/( r '~ 2 ), one 
would expect that d(T, y) is somewhat larger than N, so that c just becomes 
positive: 

c fa N « 1. 

Clearly, a lower bound for d(r, y) suffices. To compute such a bound we use 
Lemma 3.5 from [9], which we now state. 

Lemma 1. If c l5 . . . , c r _ 1 is an LLL-reduced basis for the lattice spanned 
by the columns of the matrix T, and (s%, . . . , s r _i) are the coordinates of 
y G Z r ^ 1 with respect to this basis, then 

d(T,y)>2-( r - 2 y 2 \\s r ^\\\c 1 \, 
where \\ ■ \\ denotes the distance to the nearest integer. 
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When a new upper N\ on n is found, the procedure can be repeated with 
N\ instead of N. 

As for the practical calculations, for p = 41 and p = 5581 we use K = 
10 39 , and for p = 271 we use K = 10 41 . For p = 41 the conjugates k) = 
(1,3,4) turned out to give the best results, and for p = 5581 we took the 
conjugates (i,j,k) = (2,3,4) in the case 7 = 71, and (i,j,k) = (1,3,4) in 
the case 7 = 72- For p = 271 we took the conjugate i = 2. The values of 
the entries of T and y are given in the appendix of the extended version of 
this paper [5]. 

As a result of our computations we found: 

for p = 41 : [c x | = 1.148... • 10 30 , 

d(T,y) > 1.017 
d{T,y) > 3.286 



for 5 = a: 
for 5 = a + 1: 
We infer n < N\ 



for p = 5581, 7 = 71: 



for 6 = a: 
for 5 = a + 1: 
We infer n < N\ 



for p = 5581, 7 = 72: |c 



for 6 = a: 
for 5 = a + 1: 
We infer n < N\ 

for p = 271: 1^1 : 
for 5 = a: 
for <5 = a — 1: 

We infer n < N\ 



1.148... 
|s 4 || = 0.2505. 

Ml 

= 59. 



0.0809 . 



_ Cl| = 1.123.. 
s 4 || = 0.4489..., 



S4|| 

= 23. 



s 4|| 
s 4|| 

= 23. 



0.3512..., 

^ I = 6.875.. 
0.3849..., 
0.4225..., 



2.826... • 10 20 , 
02 1| = 0.2302.. 
S2II 
= 37. 



0.2565.. 



■ 10 30 , 
d(T,y 
d(V,y 

■ 10 29 , 
d(T,y 
d(T,y 



d(T,y 
d(T,y 



10 29 , 
10 28 , 



0.0650..., 
0.0125.... 



> 1.784 

> 1.395 



> 9.357 

> 1.027 



> 4.602 

> 5.127 



10 29 , 
10 29 , 



10 28 , 
10 29 , 



10 19 , 
10 19 , 



0.1119..., 
0.0867.... 



0.0568..., 
0.0628.... 



0.0014..., 
0.0050.... 



All reduced upper bounds are well below the lower bound n > 239 we 
had already found in Section [2j Hence, the given equations have no positive 
integer solutions (n,x). 

We used the built-in LLL implementation of Mathematica 7.0. The total 
computation time was about 0.5 second on a standard laptop. 
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